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In this paper we examine gauge coupling unification of the the non-SUSY SO(IO) grand unified 
theory proposed by Babu and Mohapatra at the two loop level [I]. This theory breaks down to 
the standard model in a single step and has the distinguishing feature of TeV non-standard model 
scalars. This leads to a plethora of interesting new physics at the TeV scale and the discovery 
of new particles at the LHC. This model gives rise to testable proton decay, neutron-antineutron 
oscillations, provides a mechanism for baryogenesis, and contains potential dark matter candidates. 
In this paper, we compute the two loop beta function and show that this model unifies to two loop 
order around 10^® GeV. We then compute the proton lifetime, taking into account threshold effects 
and show that these effects place it above the Super Kamiokande limit [2]. 


I. INTRODUCTION 

It has been known for quite some time that the stan¬ 
dard model, despite being the most predictively success¬ 
ful theory in physics, is far from complete. A flagrant sign 
of this deficiency came with the observation of neutrino 
oscillation, indicating that neutrinos have mass in con¬ 
tradiction with the standard model [3HT4] . Grand unified 
theories (GUTs) pose a solution to this problem that is 
both simple and elegant. GUTs also suggest resolutions 
to many other issues, including the origin of charge quan¬ 
tization and the relative size of the standard model gauge 
couplings. In addition, they also suggest resolutions to 
cosmological problems, such as providing a mechanism 
to explain baryogenesis and candidates for dark matter. 

Among this broad class of theories, the subset of mini¬ 
mal SO(IO) theories are arguably the most natural. They 
describe each generation of the standard model in an 
irreducible representation of the gauge group and have 
inherent left-right chiral symmetry. These theories also 
provide a natural explanation for the origin of neutrino 
masses and flavor oscillation through a high scale see-saw 
mechanism [3 HMS] coming from symmetry breaking 
due to a 100126-plet Higgs sector [2^ . 

A key test of grand unified theories is through proton 
decay experiments. In the standard model, both baryon 
and lepton numbers are conserved by perturbative inter¬ 
actions and hence protects protons, the lightest baryon, 
from decay. GUTs however, have additional heavy vec¬ 
tor bosons which do lead to proton decay. By nature 
of gauge symmetry breaking, these heavy gauge bosons 
have mass on the order of the unification scale which 
this proton decay mechanism to order l/M^, rendering 
protons effectively stable. 

This places a special emphasis on proton decay ex¬ 
periments, making them a powerful tool for determining 
the viability of GUT models. In experiment, the most 
commonly sought decay modes are the p —> e“'"7r*^ and 
p —>■ . For our purposes, we will focus solely on the 
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first process as it is the dominant decay channel for this 
class of models assuming that there is small flavor mix¬ 
ing [22] . By measuring this decay, the latest experiments 
have placed the lifetime of the proton on the order of 10^^ 
years mm- This places a strong constraint on the class 
of physical theories, providing a minimum requirement 
for all GUT models. 

However, there are several factors which make the the¬ 
oretical prediction of the proton decay lifetime hard to 
estimate. There are two main sources of uncertainty in 
this calculation. The first source comes from the fact 
that quark confinement is a non-perturbative effect, so 
that we must rely on numerical lattice QGD techniques 
to approximate proton decay amplitudes. The second 
source of uncertainty, comes from the fact that the low 
energy effective field theory integrates out many heavy 
fields which become relevant at higher energies with ad¬ 
ditional unknown parameters. This uncertainty will be 
taken into account in calculating the threshold correc¬ 
tions. 


II. BABU-MOHAPATRA MODEL 

In a recent paper [T|, Babu and Mohapatra pointed 
out that there exists a non-SUSY SO(IO) GUT model 
where the seesaw scale is close to the GUT scale. But the 
inherent quark-lepton unification in GUTs implies that 
AL = 2 breaking required for seesaw give rise to observ¬ 
able AB = 2 processes leading to neutron-antineutron 
oscillation. Key to this observation is the existence of a 
TeV mass color sextet scalar (A^cd^) transforming like 
(6,1,1/3) which is necessary for one loop unification. 
Unification additionally requires the existence of 2 com¬ 
plex weak triplets, which we interpret as arising from an 
additional 45 0 45'-Higgs multiplets. This collection of 
scalars lead to one loop unification of the gauge couplings 
near 10^® GeV CQ 


^ It has also been shown recently in |23l I24| that coupling 
unification can be achieved with other TeV scalar multiplets 
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This leads to several experimentally interesting fea¬ 
tures: (i) the existence of TeV scalars which have the 
potential to be discovered at the LHC in the near future, 
(ii) the prediction of proton decay with a short enough 
lifetime to be observed in the next generation of experi¬ 
ments, (iii) the existence of a color sextet based mecha¬ 
nism for neutron-antineutron oscillation which could be 
measured in the next generation of experiments |25j , (iv) 
the existence of a color sextet based mechanism for GUT 
scale baryogenesis which is unaffected by electroweak 
sphaleron processes, and (v) the inclusion of the weak 
scalar triplets which make an excellent candidate for dark 
matter due to the 45 0 45' Higgs multiplets’ lack of 
Yukawa coupling to the standard model fields. 

Since neutron-antineutron oscillation scales as M7^ 
it is important to know the precise value of the color sex¬ 
tet scalar mass [1] . Similarly, to have a reliable prediction 
for proton decay, one needs a precise value of the unifi¬ 
cation scale. To accomplish these goals, it is necessary 
to carry out a two loop analysis of the gauge coupling 
evolution which can determine the color sextet mass as 
well as the unification scales more precisely than the one 
loop analysis of [T]. It is the goal of this paper to carry 
out this program. 


III. COUPLING UNIFICATION 

The standard method for constructing non-SUSY 
SO(IO) models is by implementing intermediate symme¬ 
try groups, each of which are broken at a different energy 
level, creating a multi-step chain from SO(IO) to the stan¬ 
dard model. A detailed two loop analysis of many SO(IO) 
breaking chains were carried out in [13 [MUM 113 ■ How¬ 
ever the breaking pattern of the model studied in this 
paper was not included: 

S'Ci(lO) —^ SU{3) X SU{2) X U{1) 

mu 

We will now take a bottom-up perspective and assume 
that in addition to the standard model, there exists a 
color sextet A„cj^c which transforms as (6,1,1^), and two 
weak triplets uJi which transform as (1,3,01H The two 
loop beta function is given by [Ml IM : 


where the different field components are implicitly 
summed over. Here, C 2 (R) and S 2 {R) are the Gasimir 
element and Dynkin index of a representation R respec¬ 
tively, and F, S, G correspond to the fermionic, scalar, 
and adjoint representation of a gauge group G. When 
G is the direct product of semi-simple terms (G = 
Gi X G 2 X ...G„), the two loop term allows for mixing 
of the gauge subgroups corresponding to contributions to 
the gauge boson propagator from field multiplets trans¬ 
forming non trivially under multiple subgroups. These 
fields add additional terms to the (3 function: 

- + E & [2S2{F,)G2{F,) + 4S2{S.)G2{S,)] ( 2 ) 

3 ^ 

where Fi is the sub-representation of F transforming un¬ 
der the gauge group Gi. Now, writing the beta function 
as: 


P{9i) = 


9i 7 I \ ' 9i 9 k 1 
(47r)2'’' + ^ (47r)4''*'= 


the numerical coefficients are given by: 
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These match with the 1-loop calculation from [I] . 

As in [1], we will also study the additional case of a 
second standard model Higgs doublet H(l,2,l/2) which 
produces the beta function: 
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Taking into account the mass of the Higgs boson and 
the threshold corrections [Ml ESHM : 




s 


(47r) 


2G2(F) + yG2(G) 


S2{F) 


( 1 ) 


-^4\C2{G)f + 


4G2(5) + -G2(G) 


S2{S) 


through more complicated symmetry breaking chains. 

^ Here we use the normalization for the hypercharge so that 
V 2 


C2{R) = [iY 


1 = ^ 

auiMu) a,{Mu) 12 tt '' ^ 

the j3i can be numerically integrated using Mathematica. 
Here we have used: 

U 9 ) = S 2 {Gv) - 21S2{Gv) log — 

Mq ^ Mp 

+ ^ S 2 {S,) log ^ + 8 ^ S 2 {F,) log 

3 3 

following the calculation from [Ml |33] where the V, Sj , Fj 
are the sets of heavy vector bosons, fermions, and scalar 
bosons respectively coupled to the gauge group factor Gi. 
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FIG. 1. Running of the standard model gauge couplings with 
a color sextet A„c^c (6,1,1/3) at 1 TeV and 2 weak scalar 
triplets uji (1,3,0). 

We found that with a single standard model Higgs dou¬ 
blet, the couplings unify at 1.42 x 10^^ GeV with = 
6.16 TeV and Ma = 1 TeV, and au{Mu)-^ = 38.5. We 
also found that with two standard model Higgs doublets, 
the couplings unify at 1.06 x 10^® with = 76.8 TeV, 
Ma = 1 TeV, and au{Mu)-^ = 38.2. 

IV. PROTON DECAY 

From our calculation for the unification scale, we can 
determine the proton lifetime in this model due to the 
decay mode: p —> In the low energy effective the¬ 

ory, this decay process is dominated by dimension 6 op¬ 
erators. These operators come from integrating out the 
heavy gauge bosons and heavy Higgs particles. However, 
from [ 1 ] we know that the standard model constrains the 
Yukawa coupling to be small enough so that the cou¬ 
pling to the scalar bosons contributions are suppressed 
relative to that of the heavy gauge bosons by several or¬ 
ders of magnitude. In general there are 5 independent 
types of such operators which lead to nucleon decay. In 
the notation of Weinberg, they are denoted [55H55] : 

oilL = id aa R’^^hFl) {^Qi-ycL 
(2) 

= (9 iaa L Qj l3bL ) {^jcFiJdr ) 

^ abed — {,QiaaLQj iSbL^iQk’ycL^dL^^oijS'y^ij^kl (^) 

iota LQjj3bL^{,Qk'ycL^ldL^^oi^'y{,'^^^ij ’ 

= id aa 

where a ,/?,7 denote SU(3) color indices, i,j,k,l are SU(2) 
isospin indices, a,b,c,d are generation indices, and L and 


® It is important to note that the mass of the scalars are not fixed. 
Once the mass is fixed for one scalar, the other is fixed by the 
condition of gauge unification. The given scalar masses were 
selected to give the maximum proton decay lifetime so that the 
scalar masses are > 1 TeV. 


R refer to the chirality. Only 4 of these (1,2,4,5) are 
relevant for proton decay. Of these only: 

^ — ^1111 — idaR'R^RiiR'yL^L d^Lh'j^^CajS'y 

^ ~ 1111 ~ idaLR0L)iRjR(^R')^aj3j , , 

' ( ^ 

^ — ^1111 — (t^aLtl/ 5 L)(ifyLt-L d^LR 

^ ~ ^1111 ~ idciR'Rj 3 R)i'^^R^R)^ald^ 

lead to the proton decay process: p —>■ 7 r°e+ [51]. Here 
we used the notation: 

^abed — y^abed ^bacd) 

We will denote the first part of these operators (the parts 
with a e~^ leptonic term) Orr' where T, T' = L, R. For 
example: Olr. — ■ In order to calculate the proton 

lifetime we need to calculate the amplitude: 

(7r°e+|C>rr'b) = -(7r°e+|C>pp, |p) (11) 

where R = L and L = R. Using this, the decay width is 
given by [551 13 S] : 


F(p —>■ 7r° -h e~^) 



^CWo*(p^7rO + e+) 

i 


( 12 ) 


where the Wq are the form factors associated to the oper¬ 
ators (i.e. VFq = ( 7 r°,e+|( 5 (*)|p)) and the C* are the 
Wilson coefficients coming from the renormalization of 
gauge couplings. These coefficients are given by running 
the coupling constants from mass scale from down to 

p [5111^: 
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Here Cf” is the coupling constant in the effective La- 
grangian at the scale /Xq and hi is from the Pi function. If 
the P function changes with energy scale, such as when 
there are intermediate symmetry stages or scalar bosons 
with low mass, then one must run the Wilson coefficients 
several times starting from the unification scale (where 
the denominator starts with Q;ff) |341 155] . 
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FIG. 2. 2-loop unification of the extended Babu-Mohapatra 
models. These are labeled as {Na,Nu,), which correspond to 
the number of color sextets and weak scalar triplets. 

For this model only the Olr and Orl operators 
are in the effective Lagrangian after integrating out the 
heavy gauge bosons. Using 2 GeV as the energy scale at 
which the proton decays (the standard in Lattice QCD 
calculations), we find that for this model: 


Cl =y4i-^ = 11.58 X 
C2=A2^ = 12.97 X 

Mu 


9u 

Mh 

ah 

Ml 


(14) 


Now that the Wilson coefficients have been calculated, 
all that is required to determine the proton lifetime are 
the QGD form factors. 

As explained before, these form factors must be com¬ 
puted using lattice QCD simulations. We will use the 
results from [37] which calculates the value: 

= {AOrl\p) = -0.103(23)(34) (15) 


where the first error is statistical and the second is sys¬ 
tematic. We now have the final formula for the proton 
decay width: 


r(p —)■ n^+e'^) = 


X 


Wp 


(Ai 




Ml 


2 


(16) 


Using the values given above for the QCD factor and the 
renormalized Wilson coefficients, we find: 

t{p ^ 7r° -b e+) ~1.1 X 10^^-o-4 yrs. (17) 


where the error comes from the QCD form factor. In 
addition to this, we will also have errors coming from 
threshold effects. A short calculation shows that thresh¬ 
old effects lead to uncertainty in the unification scale: 


A log 


M^ 

Mu 


= 0.3747?i26 + 0.01l77io 
-b0.006?745 + 0.037ry^^ 


(18) 


where rji = log . Following the prescription of |5D1 - 
|32l [38] we allow the Mi/Mu to vary between 10 and 
10“^ to give an upper bound on the estimate: 

Mu/Ml = 10±°'43 ( 19 ) 

where M^ is the central value of the unification scale. 
This leads to the proton lifetime of the order: 

t{p —>• 7r° -b e^) ~1.1 X lO^^-o yrs. (20) 

Taking into account the uncertainty from the lattice 
QCD calculation and threshold effects, we can see that 
the lifetime of the proton is outside the lower limit set 
by the Super Kamiokande experiment (1.4 x lO^"^ years 
Eiin]), but easily within reach of future experiments. 


V. EXTENSION OF THE BABU-MOHAPATRA 
MODEL 


It is interesting to note that the proton lifetime of 
this model could be lengthened by extending the Babu- 
Mohapatra model. This can easily be achieved by adding 
more non-standard model scalars (A^c^^c, uj) whose num¬ 
bers we will denote (Aa, A^) for the number of sextets 
and triplets respectively. This causes the unification scale 
to increase while coupling at unification decreases. Some 
values are plotted in Figure 2 with corresponding unifi¬ 
cation values in Table 1. 


{Na,N^) 

rrit^ (TeV) 

t{p 

-fi> 

7r° -b e+) yrs. 

(1,2) 

6.2 

1.1 

X 

1033lo:4±i-72 

(2,3) 

21 

bo 

X 

10377o:4±3-21 

(3,4) 

48 

bo 

X 

104U^1±4.72 

(4,5) 

100 

5.4 

X 

1043lo:4±6-22 

(5,6) 

220 

bo 

X 

10537o:4±7.73 


TABLE I. Comparison of extensions to the Babu-Mohapatra 
model by adding additional color sextets (A) and weak 
triplets (ca). In these models rtiA = 1 TeV. 


For these models we have fixed the sextet scalar mass 
m\ = 1 TeV in order to showcase the extensions with 
new TeV level physics. It seems that in general, full two 
loop unification with low TeV scalar masses can only be 
achieved if there is 1 more triplet than sextet (although 
it is believed that all extensions will unify with higher 
scalar masses). These extended Babu-Mohapatra mod¬ 
els provide a large parameter space that can be compared 
with experimental bounds on dark matter from experi¬ 
ments such as LUX and the LHC while maintaining a 
sufficiently long proton lifetime. Further study of these 
models is required. 
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VI. CONCLUSION 

We have demonstrated in this paper, that the non- 
SUSY SO(IO) GUT model put forth by Babu and Mo- 
hapatra in [T] unifies to two loop order at Mu = 1.42 x 
]^q15±o. 43 QgY -v^rith scalars at the TeV scale [T]. We have 
presented the calculations to show that this model pre¬ 
dicts the proton lifetime to be 1.1 x yrs, 

which places the proton lifetime above of the Super- 
Kamiokande limit [5] . In addition we have also proposed 
a class of extensions of the Babu-Mohapatra model which 
incorporate additional scalar fields and have the addi¬ 
tional merits of a longer proton lifetime and dark matter 
candidates with higher mass. These models have been 
shown to meet the minimum requirement for proton life¬ 
time and therefore should be further investigated for their 
rich low energy phenomenology. Because of the inclu¬ 
sion of low mass scalar sextets and triplets, this model 
will potentially produce many physical effects which will 


be observable in the next generation of experiments. In 
addition to the prediction of particle discoveries at the 
LHG, this model is also of great interest due to its pre¬ 
dictions of neutron-antineutron oscillation, new mecha¬ 
nism for scalar mediated baryogenesis, and dark matter 
candidates [1]. 
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